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One-electron tunneling through a quantum dot with a strong magnetic field in the direction of 
the current is studied. The linear magneto-conductance is computed for a model parabolic dot with 
seven electrons in the intermediate states and for different values of the magnetic field. It is shown 
that the dot density of states at low excitation energies can be extracted from a precise measurement 
of the conductance at the upper edge of the Coulomb blockade diamond. We parametrized the 
density of states with a single "temperature" parameter (in the so called "constant temperature 
approximation"), and found that this parameter depends very weakly on the magnetic field. 
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I. INTRODUCTION 

The experimental study of vertical transport through a semiconductor quantum dotlll have stressed the similarities 
of these small electronic systems with real atoms. Shell effects (in a confinement potential which is approximately 
parabolic in shape) and spin effects are clearly distinguished in the positions of conductance peaksH. The inclusion 
of a relatively strong magnetic field, B, jpakes possible to measure spin reordering with B up to the formation of 
a completely polarized electronic dropletB. The technique has also .sevealed its power as a spectroscopic tool in the 
determinatiou-of the low-lying energy levels of the few-electron dot^. More recent developments include integer spin 
Kondo effectsO, and the study of two-electron tunneling in the Coulomb blockade regima3. 

The position of conductance peaks, obtained from the experiments, determine addition energies and even excitation 
energies to the first excited states, if they are well separated from the rest of the states. This only happens for the 
few-electron dots and for the very first excited states. In the six-electron dot, for example, at excitation energies 
around 1 meV, the density of states may be as high as 80 levels/meV. That is, mean level distance around 0.012 meV 
(see Section III). In the present paper, we show that for relatively small dots and excitation energies < 1 meV, the 
density of states can be obtained from a precise measurement of the height and position of the conductance peak at 
the upper edges of the Coulomb blockade diamonds. 

We present model calculations for a 6-electron dot in a magnetic field 8.75 < B < 12 T. Our model parameters 
are chosen to approximate the experimental conditions in papers |^,^. The relatively strong magnetic field guarantees 
that only spin-polarized states are relevant in tunneling processes. On the other hand, the temperature is low enough 
(a few mK) for thermal excitation to be neglected. Thus a pure quantum mechanical description is used. 

The plan of the paper is as follows. The details of the model quantum dot are specified in Section II. In the 
next section, the low-lying energy levels of the 6- and 7-electron dots are computed, and the density of levels at low 
excitation energies is parametrized with a single "temperature" parameter. We show that this parameter depends 
very weakly on the magnetic field. The calculation of the linear magneto-conductance (transmission coefficient) is 
sketched in Section IV. In Section V, we show how the conductance at maxima can be related to the density of energy 
levels, and we compare the estimate following from the conductance with the actual (calculated) density of states. 

II. THE MODEL QUANTUM DOT 

A schematic representation of the vertical profile (z axis) of the bottom of the conduction band for the model 
structure to be used is given in Fig. la. This is a symmetric AlGaAs(7 nm)-InGaAs(12 nm)-AlGaAs(7 nm) quantum 
well, with n-doped source (S) and drain (D) contactsO. The quantum dot is formed withiii the well region. For the 
smallest dots, the lateral {xy) confinement is approximately parabolic, with hujo = 3 meVu. 

The source potential will be taken as the reference potential. It will be fixed to zero. The barrier height due to 
the Al concentration is estimated in Ref. as 50 meV. We will add 42 meV corresponding to the Coulomb barrier 
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(6 electrons times 7 meV). Thus, the top of the left barrier is at 92 meV, and for the right barrier at {92-Vsd) nieV 
{Vg and Vsd are in fact energies, not voltages). The drain potential is —Vsd- Notice that our Vg parameter is not the 
one used in experiments. In our model, the depth of the potential well is — (Vg + Vsd/2), while in experiments Vg is 
related to the depth of the well through the capacitance of the system. 

The electron mass and dielectric constant are taken as 0.067mo and 12.5 respectively all over the structure. In the 
computation of the A^-electron states in the dot, a first quantum well sub-band approximation is used. The energy 
of the first quantum well state, Ez, is computed numerically as a function of Vg and Vsd- The total energy is then 
written 

E = NEz+E,y, (1) 

where E^y is obtained from a two-dimensional calculation which includes the lateral confinement and the effective 
Coulomb potential (equal to 0.8 times the two-dimensional potential to account for averaging in the z direction), p 

Magnetic field values around 10 T are considered. In this region, the electronic droplet is completely spin polarizedH. 
Moreover, the separation between Landau levels (LL) is fiuic — 1.728 B mcV, where B is measured in Teslas. Thus, 
the second LL is around 17 meV above the first. If|-we are interested in excitation energies below ImeV, we may 
neglect contributions from the second and higher LL'sQ. In section III, a large basis of Slater first LL functions is used 
to construct the two-dimensional matrix hamiltonian, which is further diagonalized by means of a Lanczos algorithm. 

In section IV, we will compute the transmission coefficient of the model shown in Fig. ^ in which the potential is 
sectionally constant, unlike a real structure. 



III. DENSITY OF LEVELS AT LOW EXCITATION ENERGIES 



In the present section, we show results for the spin-polarized energy levels of the 6- and 7-electron dots. The starting 
point is the "first LL" H^y hamiltonian 

H,y = i\M\ + N)hn - \M\hu;c/2 + f3 ^ < i,j\l/r\kj > a+a+aiak, (2) 

i,j,k,l 

where ujc = eB/{mc) is the cyclotronic frequency, and 17 = ^y ujq + {hujc/2)^. The sum runs over indexes, I, rep- 
resenting angular momentum projection onto the z axis, a/" creates one electron in a harmonic oscillator state of 
frequency with zero radial quantum number and / < 0. The Zeeman energy is not included. The constant /3 equals 
0.8 e'^ /{kIq), where In is the harmonic oscillator length, and k - the dielectric constant. 

The hamiltonian H^y is diagonalized in a basis of Slater determinants with fixed angular momentum projection 
M = X^tLi^*- F'o'^ ^ given M, this basis is finite. In a 7-electron system, for example, the sector with M = —80 
contains 40340 functions. These large matrices are better diagonalized with a Lanczos algorithm. 

We show in Fig. |^ the lowest energy levels (excitation energy < 1 meV) for the N — 6 and 7-electron dots in 
magnetic fields B — 8.75 and 12 T. A few remarkable facts are evident from this figure. The average number of levels, 
n, with excitation energies < AE may be very well fitted by a "constant temperature approximation"Ll 

n = exp(AS/e). (3) 

The temperature parameter, Q, exhibits a weak dependence on the magnetic field. When B varies from 8.75 to 
12 T, for example, 1/0 for the 6-electron dot changes only from 2.99 to 3.20 (meV)~^. In the language of "filling 
factors", 1/ ~ Mo/Mgs (where Mq = —N{N — l)/2 is the momentum corresponding to filling factor one), one has 
jy = 15/35 = 3/7 a.t B = 8.75 T, and ly = 15/45 = 1/3 for B = 12 T. The gap to the first excited state shrinks to 
zero at the B values where the ground-state momentum, Mgs, changes first from -35 to -39, and then from -39 to -45. 
O, however, varies very weakly. It means that Q is not a measure of this gap, but of the actual low-lying density of 
levels. 

The average density of levels following from Eq. (^) is 

^ = lexp(A£;/e), (4) 
which may be taken as a smooth version of the actual density 

^^J2^(^^-^^r), (5) 

r 

where AE^ denotes the excitation energy of the r-th level. For = 6 and B = 12 T, for example, the average density 
is around 80 levels/meV at AE=1 meV. 
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IV. LINEAR MAGNETO-CONDUCTANCE 



We willgompute the conductance, dl/dVsd, from the transmission coefficient, T, by means of a simpHfied Landauer 
expressionH 

Only one spin polarization is considered in (j^) due to the quenching of the spin-down states by the magnetic field. T 
is computed from the usual relation between transmitted and incident current in the z direction. 

The fact that there is only one electron tunneling through the structure allows one to write an ansatz for the total 
wave function in which N electrons are permanently confined inside the dot. We restrict our attention to 6 confined, 
and a seventh tunneling electron. Thus, we shall work in the interval of Vg values corresponding to the N=6 Coulomb 
blockade diamond of Fig. ||b. We shall focus on the first conductance peak when Vsd is varied, i. e. the edge of the 
diamond. 

As mentioned above, only spin polarized states are considered. The reason is the following. Initial states for the 
tunneling processes are states with TV electrons in the dot and one free electron in the source electrode. The energy 
is given in Eq. (^). Contributions to the trasmission coefficient come from intermediate states of iV + 1 electrons 
in the dot, and final states with N electrons in the dot and one electron in the drain, which energy is less or equal 
than the energy given in Eq. (j^). As we are interested in the first conductance peak when Vsd is increased, only one 
intermediate state contributes, i. e. the ground state of A'' + 1 electrons in the dot. In a strong magnetic-field, this 
is a spin polarized state. The long spin-relaxation times in quantum dots due to phase space reductiontiliil, make 
spin-ffip processes in the ffiial stage of tunneling impossible. Thus, relevant final states are also spin polarized. 

The total wave function is written in a separable way: ^ — xy The ansatz for the z function is the following: 

6 

-^.^X{z7)\{xi{zu). (7) 

u=l 

where xi is the ffi'st quantum well function. We will determine the combination xi^i)^ xy in each z interval where 
the z-potential is sectionally constant. For '^xy, a Fock representation will be used. We will simplify the notation 
further, writing \a) instead of '^xy{N)^ and I7) instead of "^xyiN + 1). |q;o) will denote the ground state. Then, the 
ansatz for the wave function is the following 
(1) z < 0: 

5*1 — a]^\aa) {exp ikiz + bi exp —ikiz} , (8) 



where ki — y'2mei/h, and the total energy is written (apart from trivial constants) as 

E = Eo., + ^ + e,. (9) 

£1 is the initial kinetic energy of the tunneling electron. We will fix it at a small value, ei = 0.01 meV. aj is the 
(xy) creation operator for an electron with angular momentum projection Iq in S. Due to the assumed cylindrical 
symmetry, the total angular momentum, M ~ + is a conserved quantity. 

(2) < z < Lb: 

^2 = ajjao) {a2 exp -k2Z + 62 expk2z} , (10) 

where £2 = ei, and k2 — \j2m{yB — £2)/^- Vb — 92 meV is the barrier height. 

(3) Lb < z < Lb + L: 

*3 = 51 It) {"3 exp ik'^z + bj exp -ikjz} , (11) 



where kj = ,/2m{Vg + Vsd/^ + el)/h, and = ei + Ea„ + Tiujc/'^ ~ E^. 
(4) Lb + L<z<2Lb + L: 



^4 = XI I") ^'^P ^ + ^4 '^^P ^4 ' 
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where fcj = y/2m{VB ~ Vsd - e?)/?i, and ef = ei + Ea„ - E^. Finally, 
(5) 2Lb + L< z: 

*5=E"Ll")«5e'^P*^5^, (13) 

a 

where = , and fcg — y/2'm{Vsd + e")/^- The sum in ([l3| ) runs over open final state channels, i. e. 

< ei + + l^d. (14) 

In the above formulae. Lb — 7 nm and L = 12 nm are the barrier and well widths respectively. Notice that our 
ansatz for ^' respects a weak version of the Pauli principle, but it is not completely antisymmetrized with respect to 
the seventh electron. 

The continuity of the current leads to relations among the coefficients a and b in the different regions. With our 
ansatz, however, we found impossible to satisfy the continuity relations for general linear combinations of intermediate 
or final states. It means that we shall compute the transmission coefficient for each pair (7, a) independently. We 
will refer to the pair (7,0?) as a tunneling channel. We will be particularly interested in the (70,0;) channels, where 
the intermediate state is the ground state of + 1 electrons in the dot. The upper edge of the = 6 Coulomb 
blockade diamond in Fig. |l]b is related to these channels. Notice that the absence of interference between tunneling 
channels should be further reinforced by temperature effects. Note also that overlapping coefRcients of transversal 
{x,y) functions cancel out in ratios b/a, and thus in T. 

The partial transmission coefhcient, r„ is defined as usual 

T„ = lDia)/Is, (15) 

where Is is the incident current, and Id{oi) is the transmitted current when only the a channel is considered open. 
The total transmission coefficient is obtained from charge conservation arguments. If there are a few open channels, 
then IsTi is the current flowing through channel 1, Is{^ — Ti)T2 is the current through channel 2, etc. The total 
coefficient is then 

T = i- Y[ (i-r„). (16) 

open channels 



V. RESULTS AND DISCUSSION 

We show in Fig. ^a the partial transmission coefhcient, T2 corresponding to the excited state a = 2 at excitation 
energy AE2 = 0.425 meV in a magnetic field B = 12 T. This channel is closed for Vsd < 0.42 meV, as expected. 
Three situations are depicted: below resonance {Vg = 39.89 meV), maximum resonance {Vg = 39.87 meV), and above 
resonance (Vg = 39.6 meV). The asymptotic shape of the curve below resonance is typical. In Nuclear Physics context, 
it is interpreted as interference between resonance and potential scatteringlj. The maximum resonance occurs at a 
Vsd slightly higher than AE2- Note that, as the incident electron energy is fixed, the maximum of T2 is not one, but 
a value around 0.85. Above resonance, the a — 2 channel remains open but the transmission coefficient diminishes. 

The computation of T from Eq. ( p^ show results like the one drawn in Fig. where all of the open channels 
(7, a) at a given Vsd are included. Vg is 39.89 meV, and i3 = 12 T. The distance between the two peaks, 0.6 meV, 
is approximately twice the excitation energy to the first excited state of + 1 electrons in the dot. Below, we will 
focus only on the first peak. 

The first two (ground and first excited) a states contribute to this peak. Both arc above resonance. The second 
and higher excited states are below resonance. Their contribution to T is very little. Thus, there are not channels at 
maximum resonance, and the value of T at the peak is lower than the maximum ^ 0.85 for the a = 2 channel. 
With a small decrease of V^, T reaches the value TJ'^'*, meaning that there is one channel at maximum resonance. 
Notice that the distance between the first and second excited a states is only 0.12 meV. It means that only levels 
which are very close in energy may be simultaneously resonant. The width of the one-channel peak, ^ 0.05 meV, may 
serve as an estimation of the resonance interval. 

The sensitivity of the first peak maximum with the number of levels in the resonance interval may be used for an 
experimental estimation of the density of levels at low excitation energies. We prove this statement in Fig. ^ The 
number of levels in an energy interval S = 0.05 meV below the excitation energy AE is drawn along with the estimate 
(points) obtained from the conductance in the following way. 
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We follow the first peak maximum for different Vg. Its position, i. e. the value of Vsd, is identified with the 
excitation energy AE, and from the height, Tpeak, we obtain the number of levels in the resonance interval as 

n = ln(l - Tp,ak)/ ln(l - T^^), (17) 

where T^^'^ is the value of T when the first excited state a = 1 is at maximum resonance. The idea behind Eq. jl^ ) 



is that the product in Ec 



(hd) may be approximated as (1 — 



We compare in Figs. | (a) and (b) the estimation coming from Eq. (17) with the exact level distribution obtained 
before. The agreement is excellent. Abrupt variations of the density of levels are nicely reproduced as well as the 
smooth behaviour obtained from Eq. (^. This agreement proves the factibility of measuring the low-energy density 
of levels in few-electron quantum dots by means of a precise measurement of the conductance. 

In conclusion, we have computed the density of spin-polarized levels of 6- and 7-electron quantum dots at low 
excitation energies and strong magnetic fields. The density is well parametrized by a constant temperature approx- 
imation. The temperature parameter, Q, shows a weak dependence on the magnetic field. We computed also the 
dot conductance (transmission coefficient) for vertical tunneling, and showed that the conductance at the upper edge 
of the Coulomb blockade diamond is directly related to the density of levels. In this way, a procedure is suggested 
for the experimental determination of the low-lying density of states from a precise measurement of the conductance. 
Altough calculations were carried out under the simplifying assumptions of zero temperature and strong magnetic 
fields, extensions to zero or weak fields and finite temperatures are also possible. 
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FIG. 1. (a) Schematic conduction band profile for the quantum dot under study. The magnetic field is aligned with the 
current along the z axis, (b) The N=6 Coulomb blockade diamond. We compute the conductance for a fixed Vg and varying 
Vsd, as represented by the dashed line. 

FIG. 2. The low-lying energy levels of the Af = 6 and = 7 dots at B = 8.75 and 12 T. 

FIG. 3. (a) Partial transmission coefficient corresponding to the second excited (final) state at i? = 12 T. (b) Total trans- 
mission at Vg = 39.89 meV and B = 12 T. 
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FIG. 4. Number of levels in the resonance interval {AE — 5, AE) with 5 = 0.05 meV. Exact results and the estimate coming 
from Eq. (0) are compared, (a) B = 12 T. (b) B = 8.75 T. 
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